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The generation of magnetic field in astrophysical bodies, e.g., galaxies, stars, planets, is one of the 
outstanding theoretical problems of physics and astrophysics. The initial magnetic fields of galaxies 
and stars are weak, and are amplified by the turbulent motion of the plasma. The generated field gets 
saturated due to nonlinear interactions. The above process is called "dynamo" action. Qualitatively, 
the magnetic field is amplified by the stretching of the field lines due to turbulent plasma motion. A 
fraction of kinetic energy of the plasma is spent in increasing the tension of the magnetic field lines, 
which effectively enhances the magnetic field strength. Current dynamo theories of are of two types, 
kinematic and dynamic. In the kinematic theories, one studies the evolution of magnetic field under a 
prescribed velocity field. In kinematic a-dynamo, the averaged nonlinear term (u x b) (u, b are velocity 
and magnetic field fluctuations respectively) is replaced by a constant a times mean magnetic field 
B . This process, which is valid for small magnetic field fluctuations, yields linear equations that can 
be solved for a given boundary condition and external forcing fields 1 3 . In dynamic theories 4 6 , the 
modification of velocity field by the magnetic field (back reaction) is taken into account. Using a 
different approach, here we compute energy transfer rates from velocity field to magnetic field using 
field-theoretic method. The striking result of our field theoretic calculation is that there is a large 
energy transfer rate from the large-scale velocity field to the large-scale magnetic field. We claim that 
the growth of large-scale magnetic energy is primarily due to this transfer. We reached the above 
conclusion without any linear approximation like that in a-dynamo. 

There is an exchange of energy between various Fourier modes because of nonlinear interactions present in magneto- 
hydrodynamics (MHD). Since there are two vector fields u and b in MHD, the energy can be transferred from u to u, 
u to b, and b to b. Energy from a parent mode u(k) or b(k) (k represents the wavenumber of Fourier mode) is trans- 
ferred to two daughter modes with wavenumbers p and k — p. The allowed triads in MHD are (u(k), u(p), u(k — p)) 
and (u(k), b(p), b(k — p)). The net effects of all the energy transfers are constant energy fluxes from large-scale 
u to small-scale u (n^>), large-scale u to small-scale b (IIj>), large-scale b to small-scale u (II„>), large-scale b to 
small-scale b (n^< ) , and large-scale u to large-scale b (11^5 ) . The superscript and subscript of II refer to the source 
and sink respectively. All these energy fluxes are illustrated in Fig. 1. These energy fluxes are analogous to that of 
Kolmogorov's flux in fluid turbulence. Note that large-scale velocity modes are forced, as shown in the figure. 

Stanisic 7 , Dar et al. 8 , and others have given formulas for computation of the above fluxes. However, Dar et al.'s 
formalism is the most general, and they have numerically computed all the fluxes of MHD using numerical data of 
direct numerical simulation 8 . In the present paper, we compute the MHD fluxes in the inertial range using field- 
theoretic method. Our calculation is up to first order in perturbation. 

We give a brief outline of the theoretical calculation ( refer to Verma 9 for details). We write down the evolution 
equation for kinetic energy spectrum ((|u(k)| 2 ) /2) and magnetic energy ((|B(k)| 2 ) /2). We carry out our analysis 
for three space dimensions, constant mass density, and zero mean magnetic field. We also assume that (u ■ b) = 0, 
and that the large-scale velocity modes are forced. There is a net outflow/inflow of energy from a wavenumber sphere 
(say sphere of radius k ) as discussed in the earlier paragraph. These energy fluxes can be easily calculated using the 
Fourier modes. To illustrate, the energy flux from the modes inside of the u-sphere of radius k to the modes outside 
of the ^-sphere of radius fco is given by 
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where 9 stands for the imaginary part, and k = p + q. Clearly the sources of energy for the b < sphere are n^5 , , 
and . If there is a net flux of energy into large-scale magnetic energy, then magnetic energy at large-scale will 
grow, or dynamo is active. 

We calculate the energy fluxes perturbatively to first order. We assume homogeneity and isotropy for the flow. In 
the correlation functions we have included kinetic helicity (Hk = (u • oS) /2, where u> is the vorticity) and magnetic 
helicity (Hm = (a ■ b) /2, where a is the vector potential), which are defined using 
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Note that helicities break mirror symmetry. 

We focus on the fluxes in the inertial range. Based on recent theoretical 10 " 14 and numerical 15,16 evidences, we take 
Kolmogorov's spectrum for the correlation functions in the inertial range, i.e., 

C uu (k) = ^-n 2 / 3 fc- n / 3 ; C bb (k) = C uu (k)/r A ; (4) 

H K (k) = r K kC uu (k); H M (k) = r M C b \k)/k (5) 

Since both magnetic and kinetic energy spectrum are Kolmogorov-like in the inertial range, ta^k and tm can 
be treated as constants. We also use turbulent or renormalized viscosity and resistivity in our calculation. These 
quantities have been recently derived by Verma 13 ' 17 . Using the steady-state condition, we also calculate the energy 
supply from the large-scale velocity field to the large-scale magnetic field 11^5 = 11^5 + . Since the inertial-range 
energy spectrum is universal, the above energy flux is independent of the details of large-scale forcing. 

We calculate various fluxes in terms of r a,^k,^m (see Verma 9 for further details). In Fig. I we have shown 
both nonhelical (solid line) and helical (dashed line) contributions for the case {ta = 5000, tk — 0.1, tm — —0.1). 
Large ta (here ta = 5000) corresponds to a very weak magnetic field compared to the velocity field, similar to 
the early phase of galactic and stellar evolution. The choice of tk — 0.1 and Tm = —0.1 is motivated by the 
fact that both kinetic and magnetic helicities in most of the astrophysical plasmas are relatively small. The recent 
simulations of Brandenburg 18 as well as the EDQNM calculation of Pouquet et al. 4 shows that the magnetic helicity 
is negative for small wavenumbers. Since the small wavenumber contributions dominate the contributions from the 
large wavenumbers, we have taken ru < 0. 

The flux ratios shown in Fig. 1 illustrate many important results. They are: 



There is a large energy flux from large-scale velocity field to large-scale magnetic field (n^5 ) . In addition, there 

b<helical~ b<helical> 



are two other fluxes, ife. u , ,. to the large-scale magnetic field. These fluxes are responsible for the 



growth of large-scale magnetic field in the initial stage of evolution. Pouquet et al. 4 , Pouquet and Patterson 19 , 
Brandenburg 18 , and many others generally highlight n^ /ieHcai inverse transfer, and do not consider 11^5. 
Recently Brandenburg 18 ' 20 argues that large-scale magnetic energy is sustained by nonlocal inverse cascade 
from the forcing scale directly to the largest scale of the box. He relates this effect to the a-cffect. In this 
paper and Verma 9 we have computed the relative magnitudes of all three contributions for generic parameters 
discussed above, and find that all of them to be comparable, however, is somewhat higher. 

• As shown in the figure, the nonhelical component of magnetic energy flux (n£>) is forward, while the helical 
component of magnetic energy flux is negative (inverse). The overall magnetic energy flux however is positive. 

In our theoretical calculation we have assumed homogeneity and isotropy of the flow, as well as Bo = 0. These 
assumptions are likely to hold in the early stages of the galactic evolution before large structures appear. To model 
the early evolution of galaxies, we assume that the large-scales contains kinetic and magnetic energies. During this 

I£<, and helical components Il b b > helical 



unsteady period, magnetic energy(£' b ) will be amplified by the nonhelical 11^5, and helical components n b > helical and 
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We assume quasi-steady state for the galactic evolution. Our theoretical calculation performed for large ta shows 
that the energy fluxes of the right-hand-side of Eq. (^) is proportional to HE b /E u , where E u is the kinetic energy, 
and n is the total energy flux or energy supply rate. Using Kolmogorov's spectrum, we obtain 



E^^E^e.plj^^tj (7) 

where L is the large-length of the system. Clearly, the magnetic energy grows exponentially in the early periods, 
and the time-scale of growth is of the order of Lj\J E u , which is the eddy turnover time. Taking L « 10 17 fcm and 
V E u w 10km/ sec, we obtain the growth time-scale to be 10 16 sec or 3 x 10 8 years, which is in the expected range 21 . 
Hence, we are able to construct a nonlinear and dynamically consistent galactic dynamo based on the energy fluxes. 
In this model, magnetic energy grows exponentially, and the growth time-scale is reasonable 21 . 

To recapitulate, we address the dynamo problem in the light of turbulent energy flux. Our approach is very different 
from the a-dynamo picture. Our theoretical calculation, based on perturbative field theory, shows large amount of 
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energy transfer from large-scale velocity field to large-scale magnetic field. This transfer occurs in both helical and 
nonhelical MHD. We believe magnetic energy growth is primarily due to this energy transfer. Regarding magnetic 
energy flux, there is (a) nonhelical forward flux from large scales to small scales, and (b) helical inverse transfer from 
small scales to large scales. The net magnetic energy transfer, however, is positive. 

We have constructed a model of galactic dynamo based on our energy flux results. We find that magnetic energy 
grows exponentially, and our estimate of its growth time-scale is consistent with the current observational estimates 21 . 
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FIG. 1. Various energy fluxes in wavenumber space for parameter (ta = 5000, tk = 0.1, tm = —0.1). The illustrated 
wavenumber spheres contain u < and b < modes, while u > and b > are modes outside these spheres. The energy fluxes are 
from large-scale u to small-scale u (II">), large-scale u to small-scale b large-scale 6 to small-scale u (IT'>), large-scale 

b to small-scale b (n^), and large-scale u to large-scale b (11^5). The velocity fields at large-scale are forced, and the net input 
energy is 1 unit. The solid and dashed lines represent the nonhelical and helical contributions respectively. 
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